Approach to accurately measuring the speed of optical precursors 



Chuan-Feng Li,'* Zong-Quan Zhou,' Heejeong Jeong,^ and Guang-Can Guo' 
^Key Laboratory of Quantum Information, University of Science and Technology of China, CAS, Hefei, 230026, People's Republic of China 
^Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire, 03755, USA 

(Dated: January 13, 2013) 

Precursors can serve as a bound on the speed of information with dispersive medium. We propose a method 
to identify the speed of optical wavefronts using polarization-based interference in a solid-state device, which 
can bound the accuracy of the speed of wavefronts to less than 10"'' with conventional experimental conditions. 
Our proposal may have important implications for optical communications and fast information processing. 
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The advent of fast information processing in quantum net- 
work devices often raises the need for accurate control of light 
pulse propagation. The speed of the information encoded on 
the optical pulse should obey the causality This seem- 
ingly obvious fact has been experimentally proved only re- 
cently |2]. Unlike the well-defined velocity of a point particle, 
the velocity of a light pulse traveling through an optical mate- 
rial is not precisely defined. The motion of a light pulse can 
be approximated by the group velocity (vg), which is given 
by = c/{n + Lodn/dco), where c is the speed of light in 
a vacuum, n is the refractive index of the material and a> is 
the frequency of the light ISJ. It can be clearly seen that Ug 
can be greater than c and even be negative when tadn/dta is 
negative (anomalous dispersion). Experimental studies of fast 
light propagation include Refs. [2, 4-7]. To resolve the appar- 
ent contradictions between fast light propagation and the the- 
ory of relativity, optical precursors were introduced by Som- 
merfeld and Brillouin in 1914 |8-10]. The theory states that 
the front edges of an ideal step-modulated pulse propagate at 
speed c because of the finite response time of any physical 
material and no components can oveitake this wave front. The 
foreiTinners, now known as Sommeifeld-Biillouin precursors, 
are followed by the main pulse traveling at its group veloc- 
ity. This conclusion, while conceptually clear, lacks of ex- 
perimental evidence directly relating to the speed of optical 
precursors. 

Optical precursors are recently of great interest because 
they have applications in biomedical imaging flT], under- 
water communications 1 12], and the generation of high peak 
power optical pulses lll3l Il4ll . A number of theoretical and 
experimental studies have been carried out lll5l - l2Tll . Specif- 
ically, direct observation of precursors which are separated 
from the delayed main fields has been achieved with elec- 
tromagnetic ally induced transparency (EIT) in cold atoms 
ll20l l2lll . For all these works, the authors claim that the ris- 
ing edge of the input pulse propagates at the speed of c in the 
form of Sommerfeld precursors. However, considering the 
rise time of the detecting systems, a time delay of shorter than 
1 ps would not be captured if the precursors travel at a speed 
which differs little from c. 

In this paper, we propose an experimental method to pre- 
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cisely determine the speed of optical wavefronts propagat- 
ing through a dispersive medium, i.e., optically pumped 
Nd-'^:YV04 crystal by interference with a reference pulse 
traveling through a pure YVO4 crystal. If there is any speed 
difference between wavefronts in the Nd^^:YV04 crystal and 
the YVO4 crystal, the interference pattern will be significantly 
altered and reveal the influence of Nd^^ ions on the speed of 
the wavefronts. 

Nd^^:YV04 crystals (doping level, 10 ppm.) have been 
systematically inv estig ated as a candidate solid-state quan- 
tum memory ll22l - l24ll . The ''/g/a — > "^^"3/2 transition of 
Nd^^ around 879.705 nm has a narrow homogenous linewidth 
(F/, ~63 kHz) and wideband inhomogeneous broadening 
(F,„/, ~2.1 GHz) at low temperature. Nd^^ is a Kramers ion 
and thus has strong first-order Zeeman effect, which yields a 
large ground-state splitting under the application of a mag- 
netic field (10-100 GHz/T). From the inset of Fig. 1, it can 
be seen that a moderate magnetic field splits the ground states 
into two Zeeman spin levels (|1), |2)) which connect to an ex- 
cited state (|3)). Thus, this provides a A-like system ll23il24ll . 

Note that any real pulse has a step-rising front which is 
composed of the infinite spectral components, even the pure 
YVO4 crystal cannot response to it, the ideal step-rising wave 
front propagates at the speed of c independent of any medium. 
While for the light pulse produced by electro-optic modula- 
tor (EOM) with a finite rise time, the YVO4 crystal behaves 
as a wideband response medium, the wave fronts which trav- 
els at the speed of c are generally too weak for detection in 
the YVO4 crystal. In both the Nd^+:YV04 crystal and the 
pure YVO4 crystal the experimentally detectable wavefront 
should propagate at the same speed vo = c/riQ, where no is 
the refractive index for far-detuned frequency components of 
the input pulse. The refractive index ideally reaches unity 
and vq becomes c for infinite frequency components. In the 
Nd^^:YV04 crystal, the optically pumped Nd^^ ions have a 
strong dispersion of the narrow-band frequency components 
that is near the resonance line of the Nd^^^. Those components 
make up the greatly delayed or accelerated main pulse. 

The frequency of the probe light is tuned to the Nd^^ ab- 
soiption line (W31). It is a square pulse modulated using an 
EOM with a temporal width of fo = 4 x 10"^ s and a lise (fall) 
time of f,. = 0.1 X 10""^ s. The input light intensity is assumed 
to be unity in the following calculations. £(0, co) is the Fourier 
transform of the input pulse. The Nd^^:YV04 crystal under- 
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FIG. 1: (Color online. The phase (black solid line) and intensity 
(red dashed line) of light after passing through a Nd^^+: YVO4 crystal 
versus the probe detuning. The inset is an energy level diagram of 
the — > *F3/2 transition of Nd^^ under a magnetic field. 



going EIT is characterized by its linear susceptibility ||2^ 

ao 4(A + iyi2)6i 



h iaP-4(A + /ri2)(A + i5i)' 



(1) 



where oro = 4130 m ' is the on-resonance absorption coeffi- 
cient; ^0 = noW3i/c - 1.55x10^ m"'; «o is the refractive index 
of the YVO4 crystal near 879.705 nm for light polarized par- 
allel to the crystal's symmetric c axis [26]; A = w - W31 is the 
probe detuning; and is the coupling laser Rabi frequency. 
26\ is the original absorption bandwidth, which is controlled 
by the pumping procedure. Here, we choose 6\ - 2kx\ 
MHz which is larger than F;,. The dephasing rate between 
the two ground states is 712 - 2nxQ.2 kHz. These parameters 
were chosen based on experimental observations [| 23[|24ll . The 
transmitted field is given by E{zi,<jj) — E(Q, a))Exp[iki{a))zi] 
with ki((jj) - ko ^Jl +xii(^) - ko{l H-;\fi(w)/2). Via a fast 
Fourier transformation (FFT), the transmitted field in time do- 
main can be obtained. There have been some analytical solu- 
tions to this problem il9il . which fit well with the results ob- 
tained using the FFT. However, to preserve the phase of the 
field for the coming interference experiments, we use FFT to 
evaluate the integrals. Fig. 1 shows the phase of the pho- 
tons experiencing (p = Re[ki(a>)zi] - kozi with zi - 0.01 
m and the transmitted-field intensity of E xp[-2Im[ki(aj)zi]]- 
It can be seen that far-detuned light experiences a smaller 
phase shift than near-resonance light. After passing through 
the EIT medium, the rising and falling edges at f = and 
t - to should not experience group delay caused by Nd^"^ 
ions. However, the main fields are approximately delayed by 

tg^2aoZi6i/\nA^mM- 

The experimental scheme is shown in Fig. 2. The laser 
should be wavelength tunable and have a linewidth well be- 
low 1 MHz. Using a polarization beam splitter (PBS) and a 
half-wave plate (HWP), the light's polarization can be tuned 
to maximal absorption for a Nd"'^:YV04 crystal. The top 
acousto-optical modulators (AOM) shift the laser frequency 
to be on-resonance with 1^23 to serve as the coupling light. 
The middle of the AOM sweep of the laser frequency around 
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FIG. 2: (Color online) Experimental setup for measuring the speed 
of optical wavefronts. The coupling light produced by the top AOM 
should overlap with the probe light in the first Nd'+:YV04 crystal 
( top, red filled) and be damped then. The pump light for the the 
first Nd^^:YV04 crystal is produced by the middle of the AOM and 
it is controlled by a MC. The input pulse is produced by the EOM. 
The lower AOM produces the pump light for the filter. The pure 
YVO4 crystal (bottom, green filled) is used as a reference path. After 
passing through the BD, the light undergoes polaiization-dependent 
detection which is recorded on an oscilloscope (OSC). 



Wi3 serves as the pump light. The mechanical chopper (MC) 
before the BD turns the strong pump light off during the detec- 
tion cycle. The input pulse is generated using a 15 -GHz EOM 
and the light's polarization is rotated to |+) = "^(1^) + 1^)) 
by the HWP, where H and V denote horizontal and vertical 
polarization state, respectively. To show that the wavefronts 
actually do travel at the speed of vq independently of the Nd^^ 
ions, we use a reference pulse traveling through a pure YVO4 
crystal. The beam displacer (BD) separates H and V polar- 
ized light. The H polarized probe light is sent through the 
Nd^^:YV04 crystal (top, red filled) such that it overlaps with 
coupling light at a small angle. The V polarized probe light 
is sent through the pure YVO4 crystal (bottom, green filled) 
of the same length, z\- The Nd^^:YV04 crystal's c axis is 
aligned in the horizontal direction and the YVO4 crystal's c 
axis is aligned in the vertical direction to avoid birefringence. 
AH of the crystals are placed in a cryostat with a temperature 
3 K and a magnetic field 0.5 T. Then H and V polarized light 
are combined using the BD and then undergo polarization de- 
pendent detection. The MC after the BD is used to block light 
to protect the filter and detectors during the pumping proce- 
dure of the first Nd^^:YV04 crystal. The transmitted light 
intensities for the |-i-) and |-) = "^(1^) ~ 1^)) polarizations 
are recorded using an oscilloscope (OSC). The intensities are 
given by 

h(zut) = \IFT{-^[Eh(zucj) + Ev(zuco)]f (2) 
V2 

I-(zut) = \IFT{-^[EH(zuco)-Ev{zucom\ (3) 
V2 



where 



Eh(zu(^) = Vl/2£(0,w)£xp[ili(w)zi], 
Ev{z\,aj)- Vl/2£'(0, a))Exp[ikozi] 



(4) 
(5) 



and IFT means Inverse Fourier Transformation. Note that the 
phase velocity is the relevant quantity for interference, which 
is the same in a pure crystal and a low-doping crystal. If the H 
polarized field and V polarized field arrive at the same time. 
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the interference shall give an ideal near-to-zero output for /_ 
in the first nanoseconds. 

To totally eliminate the main fields with the central fre- 
quency components, which experience greater phase shifts, 
another Nd''^:YV04 crystal with length zi is placed before 
the detectors for |-) polarized light. The crystal's c axis is 
aligned in the -45° direction for maximal absorption. It acts 
as a strong absorption filter which can be characterized by 



Si 



h) -(A + idi) 



(6) 



This is the special case of Eq. (1) with Oc = 0. 62 is controlled 
by the pump light produced by the lower AOM in Fig. 2. In 
the following text, we refer to this Nd''^:YV04 crystal as the 
"filter" for simplicity. With the filter present, the recorded 
intensity for /_ is given by 



/_(z2,f) = \IFT{-^[Eh{zuu) - Ev(.Zuco)]e' 
where k2{u)) - h){\ +xi{u>)l2). 



(7) 
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FIG. 3: (color online). Intensity of the transmitted light of the po- 
larizations |+> (black solid line) and |-) (red dashed line) without the 
filter (a) and with the filter (b). 



Z\ - 0.01 m and zi - 0.01 m for panel (b). The input pulse 
is the same as before and 62 was chosen to be 1505i. For Fig. 
3(a), we first note that the V polarized light remains the same 
as the input square pulse after traveling in the pure YVO4 crys- 
tal as mentioned above. When the H polarized wavefronts 
decay, the V polarized light starts to dominate. Therefore, 
/+ (black solid line) and /_ (red dashed line) have nearly the 
same intensity. Then, the H polarized main fields slowly in- 
crease. /+ rises with it and decays, because the delayed 
main fields are still almost entirely in the phase with the V 
polarized light. After 4 //s, the V polarized pulse ends, and /+ 
and /_ share exactly the same intensities. It can be seen that 
/+ rises immediately after f = 0, but /_ is suppressed during 
the first nanosecond. This shows that H polarized wavefronts 
travel at the same speed as the V polarized light in the pure 
YVO4 crystal. If this were not the case, /_ would rise at the 
same time and with the same intensity as /+, because only H 
(V) polarized light would present and there would be no inter- 
ference between the two paths for the given delay time. A time 
delay T4 between the arrival of H and V polarized light would 
lead to a remarkable intensity for which can be estimated 
to be Tcilty. However, for Fig. 3(a) with no filter present, /_ 
can grow to 0.004 in the first nanoseconds with no delay for 
either path. Therefore, it becomes impossible to measure 
with subpicosecond accuracy. 

The transmission /_(z2,0 with the filter is shown in Fig. 
3(b). /+ is the same as that in Fig. 3(a) and is not shown in 
Fig. 3(b). Because the filter has wideband absorption, only 
photons with frequencies far-detuned from the resonance are 
transmitted. The main fields are totally absorbed by the filter, 
and only rising and falling edges at f = and t - fo remain. 
For all time scales, at most /_ can grow to 5 x 10"^ in the 
first nanoseconds after t - Q and t - t^. Because the re- 
maining far-detuned light experiences little phase shift, there 
is near-perfect interference between the two paths. Therefore, 
whenever there is an intensity larger than 10"^ recorded for 
we may infer that there is a time delay between the arrival 
of H and V polarized light. Moreover, the interference yields 
the same results when the first Nd"'^:YV04 crystal acts as a 
two-level system, which is the fast-light regime. 
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To show the advantages of having the filter present, the sim- 
ulation results using FFT with 0.01-ns sampling resolution are 
shown in Fig. 3 with zi - 0.01 m and Z2-Q for panel (a) and 



FIG. 4: Intensity of |-> polarized light for = (a), = -40 fs (b), 
Trf = 40 fs (c) and = 4 ps (d). 
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To accurately determine the speed of wavefronts using this 
method, let zi = 0.055 m, zi = 0.08 m, and 62 = 10505i 
which is roughly the inhomogeneous broadening of the Nd^+ 
ions. Thus, the lower AOM in Fig. 2 may not be neces- 
sary. The propagation time in the first crystal is ti ^ Zi/i/o - 
400 ps. A time delay Td between the H polarized light and 
V polarized light is artificially introduced. Let Enizi, w) = 
Vl/2-£'(0, cj)Exp[iki{a>)zi]Exp[ioL>Td]- This means the H po- 
larized light will fall behind (move ahead of) of the V po- 
larized light by |tj| for > (t^ < 0). The output inten- 
sities calculated using FFT with 0.1-ps sampling resolution 
with T4 = 0, -40 fs, 40 fs and 4 ps are shown in Figs. 4(a)-(d), 
respectively. Clearly /_ with |r^| = 40 fs is much greater than 
that without delay. The longer the delay, the greater the I 
that can be obtained. Whether the H polarized light is faster 
or slower than the V polarized light has little effect on the out- 
put. This is shown in Fig.s 4(b) and 4(c). The sign of may 
be experimentally identified by introducing a compensation 
delay in either path. 

Considering a realistic detector with a bandwidth of 5 GHz, 
the recorded power for |-) polarized light can be estimated 
by averaging the results in Fig. 4(b) for every 0.2 ns. The 
maximum outputs are 9.2 x 10"", 2.5 x 10^^,5.9 x 10"'" and 
6.6 X lO"*" for Figs. 4(a)-(d), respectively. The power of probe 
light can be chosen to be 1 mW in the experiment. Therefore, 
when the power of |-) polarized light is recorded with a value 
over 2.5 pW, we can conclude that there is a speed difference 
of Au > i»o/10^ for the two paths, which is approximately 
Av/vq ^ Tii/ti. With the parameters chosen here, a bound of 
less than lO"'* on the accuracy of the wavefronts' speed can 
be experimentally determined. 

The highest speed resolution with this scheme is not limited 
to 10""^ which depends on several parameters. The first is the 
rise time of the input pulse, which is determined by the modu- 
lation method. Because /_ oc Tj/tr, the smaller is, the more 
sensitive /_ is to r,;. The faster the signal rises, the harder they 
can "see" Nd-*^ in the Nd^^:YV04 crystal. The second is the 
length of the first crystal, zi, which determines the transmis- 
sion time, fi. Because Av/vq - r^z/fi, for a larger fi, the speed 



difference can be determined more precisely. Third, detectors 
with better low-intensity responses can resolve a smaller r^. 
Note that the rise time of the detectors is not so important in 
this scheme because only the maximal intensity is of interest 
rather than the absolute time that the signal rises. Therefore, 
this method can determine the speed of the precursor to an ac- 
curacy of a given lower bound with better experimental con- 
ditions. Compared with previous works which determine the 
speed of information to an accuracy of about 0. 1 by direct in- 
tensity recording I^IItII, our scheme shall give more accurate 
measurements. 

To summarize, we propose an interference scheme for ac- 
curately measuring the speed of wavefronts in an optically 
pumped Nd^^;YV04 crystal. The simulation results computed 
with linear susceptibility of FIT medium show that the arrival 
time of wavefronts is independent of the strong dispersions 
caused by Nd^*^ ions. Considering the current available mod- 
ulation and detection technologies, the pure YVO4 crystal is 
treated as an infinite fast-response medium like the vacuum, 
so the wavefronts' speed we are measuring is wq. In principle, 
with a powerful wideband filter and single-photon detectors 
with extremely low noise, the precursors which travel at the 
speed of c in a Nd:YV04 crystal also can be identified with 
our scheme. This scheme can be carried out in other systems 
where no background wideband medium is present, such as 
cold atoms and room temperature atoms. The conclusion can 
be made general as the ideal step-rising wave fronts travel at 
exactly the same speed c in all materials. This scheme pro- 
vides an approach to accurately measure the arrival time of 
precursors in various systems with realistic detecting systems. 
Future experimental measurements will provide a strict test 
of whether the information velocity with medium violates the 
Einstein causality. This work bridges the gap between the 
ideal concepts of precursors and the real pulses with finite rise 
(fall) time by showing an accurate boundary. 
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